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FOREWORD 


"Thermal radiation" is electromagnetic radiation emitted by matter ina 
state of thermal excitation. The energy density of such radiation in an en- 
closure at constant temperature io given by the well known Planck formula. 
The importance of thermal radiation in physical problems increases as the 
temperature is raised; at moderate temperatures (say, thousands of degrees 
Kelvin) its role is primarily one of transmitting energy, whereas at high 
temperatures (say, millions of degrees Kelvin) the energy density of the radi- 
ation field itself becomes important as well. If thermal radiation must be 
considered explicitly in a problem, the radiative properties of the matter 
must be known. In the simplest order of approximation, it can be assumed 
that the matter is in thermodynamic equilibrium "locally" (a condition called 
local thermodynamic equilibrium, or LTE), and all of the necessary radiative 
properties can be defined, at least in principle. Of course whenever thermal 
radiation must be considered, the medium which contains it inevitably has 
pressure and density gradients and the treatment requires the use of hydro- 
dynamics. Hydrodynamics with explicit consideration of thermal radiation is 
called ''radiation hydrodynamics". 

In the past twenty years or so, many radiation hydrodynamic problems 
involving air have been studied. In this work a great deal of effort has gone 
into calculations of the equilibrium properties of air. Both thermodynamic 
and radiative properties have been calculated. It has been generally believed 
that the basic theory is well enough understood that such calculations yield 
valid results, and the limited experimental checks which are possible seem to 
support this hypothesis. The advantage of having sets of tables which are 
entirely calculated is evident: the calculated quantities are self-consistent 
on the basis of some set of assumptions, and they can later be improved if 
calculational techniques are improved, or if better assumptions can be made. 
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The origin of thie set of books was in the desire of a number of persons 
interested in the radiation hydrodynamics of air to have a good source of 
reliable information on basic air properties. A series of books dealing with 
both theoretical and practical aspects was envisaged. As the series materalized, 
it was thought appropriate to devote the first three volumes to the equilibrium 
properties of air. They are: 


The Equilibrium Thermodynamic Properties of Air, 
by F. R. Gilmore 


The Radiative Properties of Heated Air, 
by B. H. Armstrong and R. WwW. Nicholls 


Tables of Radiative Properties of Air, 

by Lockheed Staff 
The first volume contains a set of tables along with a detailed discussion of the 
basic models and techniques used for their computation. Because of the aize of 
the related radiative tables and text, two volumes were considered necessary. 
The first contains the text, and the second the tables. It is hoped that these 
volumes will be widely useful, but because of the emphasis on very high tempera- 
tures it ia clear that they will be moat attractive to those concerned with nuclear 
weapons phenomenology, reentry vehicles, etc. 

Our understanding of kinetic phenomena, long kncwn to be important and at 
present in a state of rapid growth, is not as easy to assess as are equilibrium 
properties. Severe limitations had to be placed on choice of material. One 
volume is offered at this time: 


Excitation and Non Equilibrium Phenomena in Air, 
by Landshoff, et al. 


It provides material on the more important processes involved in the excitation 
of air, criteria for the validity of LTE and special radiative effects. 

A discussion of radiation hydrodynamics was felt to be necessary and another 
volume was planned to deal with this topic: 


Radiation Hydrodynamics of High Temperature Air, 
by Landshoff, Hillendahl, et al. 


It is not ready for publication at this time. It will review the basic theory of 
radiation hydrodynamics and discuss the application to fireballs in the atmosphere. 
The choice cf material for these last two volumes was made with an eye to 

the needs cf the principal users of the other three volumes. 
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Most of the work on which these volumes are based was supported by the 
United States Government through various agencies of the Defense Department 
and the Atomic Energy Commission. The actual preparation of the volumes 
was largely supported by the Defense Atomic Support Agency. 

We are indebted to many authors and oryanizations for assistance and we 
gratefully acknowledge their cooperation. We are particularly grateful to the 
RAND Corporation for permission to use works of F, R. Gilmore and H. L. 
Brode and to the IBM Corporation for permission to use some of the work of 
B. H. Armstrong. Most of the other authors are employed by the Lockheed 
Missiles and Space Company, in some cases as consultants. 

Finally we would like to acknowledge the key role of Dr. R. E. Meyerott 


of LMSC in all of this effort, from the initial conception to its realization. 


We are particularly grateful to him for his constant advice and encouragement. 


Criticism and constructive suggestions are invited from all readere of 
these books. We understand the: much remains to be done in this field, and 
we hope that the efforts represented by this work will be a stimulus to its de- 


velopment. 


The Editors 
J. L. Magee 


H. Aroeste 
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Preface 


This volume is concerned with the thermodynamic properties of air 
as well as individual air constituents from 1000°K to 1070, for densities 
between 10 and 107? times sea level density. It consists of text and 
tables which were both prepared by Dr. Forrest R. Gilmore. He is also 
the author of a somewhat smaller work (Gilmore, 1955) on the same topic 
which has received wide circulation, 

The text describes the techniques, models and approximations used 
in calculation of the tables. It is not intended to be a general treatment 
of the theory of such calculations which is too well known to need restate- 
ment. The reader is expected to have some familiarity with both thermo- 
dynamics and statistical mechanics, Discussion of the approximations is 
detailed enough so that the accuracy can be readily assessed and improve- 
ments can pe readily considered. 

The tables can be accepted as an accurate “ummary of the thermo- 
dynamic properties of air and its constituents. They have been compared 
with other calculations (which differ in various approximations) and, where 
possible, with experiment. Of all inputs for calculations of hydrodynamics 
and radiation hydrodynamics, the equilibrium thermodynamic properties are 
by far the most reliable. It is unlikely that there will be a significant 
improvement of our knowledge in this area in the near future, nor is it needed, 
except perhaps at the highest density considered. 

Information provided in the tables is much more widely useful than for 
air proplems alone. It should pe noted that most of the tables are devoted 
to the individual air constituents. 

The air of this volume is a mixture of about 78% No. 21% Oo. 1% Ar, 


with a trace of CO, . It is essentially the air (i.e., dry air) found in the 
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homosphere, or the atmosphere below 95 kilometers altitude. In the 
homosphere, mixing processes are rapid, and tha composition remains 
essentially conatant. The local thermodynamic equilibrium (LTE) approxi- 
mation of radiation hydrodynamics is reasonable only in the homosaphere. 
For convenience of the reader a discussion of properties of the atmosphere 
prepared by A. D. Anderson is included as Appendix A. 

We would like to thank Dr. Forrest R. Gilmore for his splendid 
cooperation and the RAND Corporation for permission to include this work 
in our series on “Thermal Radiation." Thanks are also due Mr. A. ND. 


Anderson for preparation of Appendix A. 


J. L. Magee 


H. Aroeste 
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CHAPTER 1. THE EQUILIBRIUM THERMODYNAMIC PROPERTIES 
OF HIGH-TEMPERATURE AIR: DISCUSSION 


1.1 Introduction 


In order to calculate the behavior of nuclear fireballs and of 
hypersonic missiles and meteorites, one needs values for the thermo- 
dynamic properties of air over a wide range of temperatures and pressures 
or densities. At temperatures above about 2000°K there are great experi- 
mental difficulties in measuring these properties directly, while it is 
generally accepted that careful theoretical calculations can yield results 
of high accuracy. Consequently, in this chapter the methods and data 
available for the theoretical computation of such properties are outlined, 
and published results reviewed briefly. Some improvements over existing 
treatments are also indicated. Extensive tables based on these improved 
expressions are given in Chapter 2. 

Throughout this chapter the assumption of local thermodynamic 
equilibrium (LTE) is made. In some situations of interest, such as fire- 
balls or missile trails at very high altitudes, this assumption is not valid, 
and the results given in this chapter are not applicable. In such non- 
equilibrium situations the thermodynamic properties of air depend upon 
the energy deposition mechanisms and the subsequent atomic and molecular 
processes, as discussed in a companion volume of Thermal Radiation 


Phenomena, Excitation and Non~equilibrium Phenomena in Air. 


; 
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1.2 Air as a mixture of idee] gases in chemical equilibrium 
The thermodynamic properties of most gases at low and moderate 
densities can be approximated Gver a certain temperature range by a 


thermal equation of state of the form 

FV = NRT (1. 2-1) 
and a caloric equation of state of the form 

5 = NCT, (1. 2-2) 


where p isthe pressure, V the volume, N _ the numberof moles, R 

the gas constant per mote, T the temperature, EE the internal energy, 

and C,, the (constant) molal specific heat at constant volume. A gas which 
satisfies Eq. (1.2-1) and Eq. (1.2-2) is called a perfect gas. Many gases 
satisfy Eq. (1.2-1) over a fairly wide range of temperatures, but Eq. (1.2-2) 
only over a much narrower range. Consequently, the definition of a perfect 

gas is sometimes modified to include gases which obey Eq. (1.2-1) but not 

Eq. (1.2-2). However, to avoid ambiguity it seems preferable to follow thermo- 


chemical practice, arid use the phrase ideal gas tc denote a gas which ts 


thermally perfect but has an internal energy varying arbitrarily with the temperature. ~ 


Ideal-gas thermodynamic properties fer many pure elements and 
compounds have beer: calculated and tabulated by various workers. The 


most extensive and accurate set of such tables was recently issued by 


a Joint Army-Navy-Air Force (JANAF) Thermochemical Panel in loose leaf 
form (JANAF Thermochemical Tables, 1960), with revisions issued 
periodically. These tables include values for most air molecules to 
cooo°rK, Methods for making such calculations at still higher temperatures 
will be discussed in later sections of this chapter. 

Air is a mixture of nitrogen and oxygen, with small amounts of argon, 
carbon dioxide, and rare gases, and a variable amount of water vapor. 

At low and moderate d2nsities {ts thermodynamic properties may be obtained 
simply by adding the ideal-gas contributions from its components. This is 
straighttorward at low temperatures, where the composition does not change. 
At high temperatures, however, dissociation and other chemical reactions 
cause the composition to vary with temperature and pressure. Because of 
this variation the product pV for the mixture is no longer proportional 

to T , even though each component ts effectively ideal. Consequently, 
aerodynamicists often call air at high temperatures a "real gas," although 
thermochemists prefer to reserve this term for high densities where inter- 
molecular forces cause deviations from ideal-gas behavior. 

The first step in calculating the thermodynamic properties of high- 
temperature air is usually to calculate the equilibrium chemical composition. 
The basic condition for chemical equilibrium (Epstein, 1$37) is that at a 
fixed temperature and pressure the amounts or concentrations of the various 


chemical species must be such as to minimize the Gibbg free energy of 


the mixture, subject to conservation of the chemical elements in the mixture. . 


If the temperature is high enough to produce significant ionization, electrons 
and ions must be included among the species, and charge must also be 


conserved. 


: 
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The Gibbs free energy of an ideal-gas mixture is given by 


O : 
Prot ™ Dy N, (FP + RT In) , (1.2~3) 
i 


where N, is the number of moles of the i-th species, FP 


gas molal free energy of this species (at the temperature T and one 


is the ideal- 


atmosphere pressure), and P, is its partial pressure (in atmospheres). 
White, Johnson, and Dantzig (1958) have developed numerical methods 

for minimizing Eq. (1.2-3) while keeping the temperature and the total 
Pressure fixed and conserving mass. Their methodr are convenient for 
calculating equilibrium compositions at specified temperatures and pressures. 
For fireball applications it is preferable to make computations at specitied 
temperatures and volumes (or densities), because the air density within a 
fireball varies by two or three orders of magnitude at most, while the 
pressure variation can be several more orders of magnitude. A convenient 
procedure for such computations can be based on the thermodynamic principle 


of minimizing the Helmholz free energy (or work function) 


Rioe  Fioy @ BY -» N, CFP - RT +RT in (N,RT/V)] , (1. 2-4) 
i 
at a fixed temperature and volume (Epstein, 1937). Since Eq. (1.2~-4) has 
nearly the same form as Eq. (1.2-3), the procedures of White, Johnson and 
Dantzig (1958) can be used with only minor modification, 
Fer systems involving a limited number of reactions, or for more 


complex systems if a high-speed computer ts not available, a computation 


method based on “equilibrium constants" may be more convenient. To 
derive the necessary equations from the free-energy minimization principle, 
consider the effect of small changes 6N, in the amounts of the chemical 
species, while the temperature and (total) pressure are held fixed. By 
the use of Eq. (1.2-3) and the relation p,* pN,/r N, . the variation 


in free snergy at constant Pp may be expressed as 


. - 
OF ot >> (F; + RT In p,) 6N, +>) N,RT 6(in Nj} >> N,RT san) Ny). (1-2-5) 
i i i 


Since § Inx= 8x/x , the last two terms on the right hand side of 


Eq. (1.2-5) cancel and the final result becomes 


Cc - 
SFrop 7D, (PEt RT In py) ON, (12-6) 
1 


In equilibrium the free energy has a minimum, so SF ot = 0 and 


from Eq. (1. 2-6) 


yy Inp, = - 2s 8N, (FV/RT) (1.27) 
i i 
for all sets of SN, which satisfy the mass-balance constraints. One 


way of satisfying these constraints is to consider possible chemical 


reactions, such as 


xY=X+yY. (1, 2-8) 
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A composition variation due only to this reaction has SNyy =z = bNy == SNy 


with all other 6N, vanishing. Eq. (1. 2-7) then gives 


FO. - FO - FS 
In Pyy -Inpy- Inpy = - RT (1. 2-9) 


Taking the exponential of this equation, one obtains 


Pp : 
oy Py = exp ~ Aik yt 3 (1.2-10) 


The right hand side of Eq. (1.2-10) is independent of the composition 
and pressure of the mixture (although it does vary with temperature); it is 
conventionally called the equilibrium constant, K, , Of reaction Eq. ().2-6). 
With the help of the ideal-gas relation, p, = N,RT/V = R,RT , where nN, 
is the concentration (moles per unit volume) of the species i , Eq. ().2-10) 


may be transformed to 


N FO, - FO - FO 
xy. OXY 
rs, = RT exp ]- RT ‘ (1. 2-11) 


(In Eq. (1.2-11) the Fy are conventionally evaluated at 1 atmosphere 
pressure; hence, the units of RT must be atm-volume/mole.) The right 
hand side of Eq. (1.2-11) may be called the concentration equilibrium constant 
K, , for reaction Eq. (1.2-8). 

If there are I chemical species present in a reacting mixture and J 
mass-balance conditions to be satisfied (one for each element present, plus 


one for charge neutrality if ions are present), it turns out that I-J 


independent equilibrium equations of the form Eq. (1.2-10) or Eq. (1.2-11) 


can be found, and thus there are just enough equations to determine the. 
I partial pressures or concentrations uniquely. The set of equations 
used, however, is not unique, since, for example, the sum or difference 
of two reactions is also a permissible reaction; its equilibrium constant 
ig Just the product or quotient of the constants for the two reactions. Any 
convenient complete set of equations can be used to get the equilibrium 
composition. 

Uniess only two or three reactions are involved, a closed-form 
sOlution to the equilibrium equations is usually not obtainable, and 
numerical iteration schemes must be used. An important exception occurs 
however, when only atoms, atomic ions, and electrons are involved. 


The equilibrium ionization equations for each element X_ can be written 


Nye = Kye Ny Mee Nye * Keen BBS) Nyaa = Kyte RyRy eee (162-12) 


*” 
where N. is the electron concentration. The sum of these concentrations, 
plus Ny , is the total concentration of the element X , Riot x + which 


is usually specified. The summed equations can be readily solved for Ny: 


N 


N = ot. 


x 2 e (1.2-13 
1+ Ky,/N, + Ky NS t eee 


If the negative ion X is stable, an additional equation for its 
concentration must be included in Eq. (1.2-12), and corresponding 
additional terms added to subsequent equations, but this does not 
complicate the analysis significantly. 


Hthe elements X,Y, ... are present, the charge-balance relation is 


Rt Ryu * Ray ty te + My, + UR, + Ny tee 6 (12-14) 


With the help of Eqs. (1.2-12), (1.2-13), and similar relations for Y , ; 
etc., Bq. (1.2-14) can be written 


t 


n Ree Ke + 2hy/Ne + ees) % Rot x he, + 2K44/Re + ore) 
© 14 KG /N, + Ry /Me +... 1+ KN, + Ky /Re +... Sites 


In @ specified-density problem the total concentration of each element is 


known, and Bq. (1.2-15) can in principle be solved for N, , although 
analytic solutions are possible only when the number of different ions present 
ig small. However, if the elemental composition (hence, the ratios 
Riot.v/ Ntot.x ' tee) and the electron concentration N @ are taken as 
independent variables, Eq. (1.2-15) can be solved directly for Niot.x ’ 
and then the other concentrations can be obtained from Eq. (1.2-12), atc., 
all without iteration.” 
Although a solution of Eq. (1].2-15) at a specified density (i.e., specified 
ulues of Nrot.X ; Biot. » @tc.) usually requires an iteration process, it 
@ particularly simple one, involving only the variable XN, « Anumerical 
study by the writer shows that a scheme using the total atom concentration 


as the first approximation to RN. + and obtaining successively better values 


This fact has been discovered, apparently independently, ody a number 
of different workers. The earliest publication known to the writer is 
that of Hilsenrath, Green, and Beckett (1959). 


| 
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by substituting the previous value in the right hand side of Eq. (1.2~15), 
converges quite rapidly when the second approximate N, is larger | 
than about half the atom concentration. For smaller N, values, 
rapid convergence can be obtained by using for each new approximation J 
the geometric mean between the previous approximation and the value | 
obtained from Eq. (1.2-15). 

Once the equilibrium composition of a reacting gas mixture such 
as high-temperature air is obtained by the free-energy-minimization or 
equilibrium=constant methods described above, the basic thermodynamic 
properties may be obtained simply by adding the ideal-gas contributions 
from the various constituents. The ideal-gas entropies and free energies 
ate conventionally tabulated for 1 atmosphere pressure, and must be 
corrected to the actual partial pressure of each species, as earlier ‘ 
indicated for the free energy (see Eq. (1.2-3)). 

Calculations for “differential” thermodynamic functions, such 
as the specific heat, are less straightforward. A specific heat computed 
by averaging the specific heats of the component species, weighted 
according to their equilibrium concentrations at a given temperature and 
density, is a "frozen composition" specific heat, useful only in situations 
where the temperature is varied so rapidly that chemical reactions do not 
have time to take place. In the usual true-equilibrium situation, however, 
the variation of composition with temperature gives an additional contribu- 


tion to the specific heat. Hochstim (1962) has presented rather lengthy 
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equations for the specific heat and related quantities in terms of the 


equilibrium constants and their temperature derivatives, and shown how 


JS ooo ~ lipase rote o 


they can be solved for air using a high-speed computer. In most cases, 
however, it is probably simpler to calculate the energy at two or three 
closely-spaced points and obtain the specific heat by numerical 
differentiation. Moreover, in many fireball calculations the energy 

and pressure of air are approximated by analytic functions of temperature 
and density. Differentiation of these functions gives approximate 
values for the specific heat, velocity of sound, etc., rather 

simply. 

The equilibrium compositions obtained in the course of calculating 
thermodynamic properties are also useful in determining radiation properties 
(see Volume 2) and transport properties (viscosity, electrical conductivity, 
etc.) of air. For these purposes it is often more convenient to express 
species concentrations in terms of molecules (atoms, ions) per unit volume. 
Since thermochemical tables are usually based on the gram mole, the 
equations in this chapter are given in molal form, but they are equally valid 
if concentrations are expressed in molecules (or particles) per unit volume, 
while volumes, energies, etc., are taken per molecule instead of per 
mole and the gas constant per moie, R , is replaced by the gas constant 


Per molecule, k (Boltzmann's constant). 


1.3. Ideal-gas properties for monatomic gases 

As will be shown below, the thermodynamic properties of monatomic 
gases depend in part on the quantum-mechanical properties of their 
constituent atoms. It is beyond the scope of this chapter to explain the 


quantum mechanics of atomic structure (for a good introduction, see 
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Herzberg, 1944). However, in the next few paragraphs, enough of the 

terminology will be explained to permit the reader to use existing energy- 

level tables to calculate thermodynamic properties. | 
An atom or atomic ion may exist in many different states with different 

energies, corresponding to different arrangements of the orbital electrons. 

According to quantum-mechanical principles for light atoms (Russell-Saunders coupling, 

Herzberg, 1944), the states are grouped in terma, where for eech term the 

total orbital angular momentum (in atomic units) of the electrons must be 

an integer. The terms are designated, for histerical reasons, by the 

letters S,P,D,F,G,H, eee , corresponding to total orbital angular 

momenta, L, equalto0,1,2,3,4,5, e+ , respectively. In 

addition, each electron possesses a spin (intrinsic angular momentum) of 


1/2 , and these spins add algebraically to give a total spin §. The 


quantity 28+ 1 is called the multiplicity of the term; it is conventionally 
written as a superscript. Thus aterm with L=2 and S=3/2 ts 


4 
4p - Agiven atom or ion may have many terms of the same 


denoted by 
type, with the same L and S§ values, but differing in the arrangement 
of the individual electrona, and thus in the energy. 

Aterm having S=s<L is composedof 28+1 Jevels, having total 
angular.momenta J# L+5S, L+S-1, L+8- 2, «e+, L- &, respectively, 
all with energies fairly close to each other. (If Ls 8S , the quantities L ! 
and S_ should be interchanged in the preceding statement). Thus, a 4p : 
state is composed of four levels designated by subscripting their J values: 
“Dap P “Divs ’ “Da /2 ‘ ‘Di, . By placing the atom in a magnetic 
field, each of these levels can be split further, into 2J+ 1 individual 
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subdlevels. Quantum theory shows that there are just this many distinct 
“quantum states” in such an electronic level; in the absence of a magnetic 
or electric field their energies are the same, 80 that they can be treated 
as a single level with a “degeneracy” or atatistionl weight of 2J+1. 
Often, the splitting of an L,S_ state into levels with different J 
values can also be ignorec, and it can be treated as a single level 
with a statistical weight of (27.+ 1) (28+ 1). Thus, a “D state 
has a statistical weight of 20; its four levels (4p, /2' ‘p, /2! ‘p, /2' 
‘D, iy have statistical weights of 8,6, 4, and 2 , respectively. 
Except for one-electron atoms and ions, accurate calculationa of 
energy levels are difficult, while very accurate measurements by 
spectrography are usually relatively simple. Since spectroscopists 
measure wavelengths, which by Planck's law are inversely proportional 
to the energy level differences, they conventionally present energy levels 
in units of om7?, To convert to true energy units, these must be multiplied 
by hc , where h_ is Planck's constant and c 1s the velocity of light. 
Moore (1949) has collected together the best experimental energy- 
level values for the light elements and ions; her tables ere complete enough 
for many practical purposes. If necessary, Moore's values can be tested 
for completeness by comparison with the tables of predicted terms given 
in the front of the same volume, and missing values estimated by using 
various semiempirical methods (isoelectronic extrapolation, quantum 
defects, screening constants (Herzberg, 1944; Edlen, 1964)). For the 
highly tonized species, for which experimental data are incomplete, Edlen 


(1964) gives a number of useful extrapolation formulas and tables. 
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Fig. 1.1 depicts, as an example, the energy levels of the nitrogen 
atom, Since, there are many terms withthe same L and 8 values, 
the terms are also marked with the quantum numbers of the individual 
electrons. Thus, 2s 2p4 indicates one electron with principal quantum 
number n= 2. and orbital angular momentum 4#0 , and four electrons 
with n= 2 and 41. Although not marked, all of the states shown 
have two inner, tightly bound electrons, is? ; states without these 
electrons lie so high in energy that they can be completely ignored for 
present purposes. 

According to statistical mechanics (Mayer and Mayer, 1940), in an 
ideal gas in equilibrium, the fraction of an atomic or ionic species which 


is in a particular electronic state or level j is 
N,/N= (9;/Q) exp (-hcB*,/kT) ‘ (1.3: 1) 


where 9; is the statistical weight of the level, E'; is its energy 

(in om” 1), and Q_ is a proportionality constant called the partition 
function. The fraction N,/N is also called the "fractional population" 
or “occupation number." When Eq. (1.3-1) is summed over all j values, 
the left.hand side reduces to unity and the equation can be solved for the 


Partition function: 


Qs 2s Qy exP (-c,E'//T) ' (1-3-2) 
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where Cy is known as the secord radiation constant, 
Cy ™he/k = 1,43879 om-°K (1.3-3) 


The values quoted in this chapter for the pertinent physical constants 
are those recently recommended by the National Academy of Sciences - 
National Research Council (1964); they differ slightly from thosa used in 
previous work. 

The electronic~energy contribution of the atoms in state j to 
the internal energy of the gas is N, G hcE' 4! where ( is the number 
of atoms per mole (Avogadro's number). Since G = R/k, this contribution 
can also be written N Re, E 1° The ideal-gas internal energy per mole 
can then be obtained by using Eq. (1.3-1) to evaluate N; ,» summing the 
contributions over J , and adding (3/2)RT for the contribution from the 


translational motion of the atoms. The result may be written in dimensionless 


form: 


are = 4 + a >> 2k" y exp (-c,E' /) ; (\.3-4) 
j 


where the superscript o indicates the ideal-gas state, the subscript 


0 indicates O°K , and 


R= 8.3143 x 10” erg/K-mole = 1.98717 cal/K-mole . (1.3-5) 


i4 


& in Eq. (1.3-4) represents the internal energy of the ideal gas at 
0°K ; its numerical value is somewhat arbitrary, since in thermodynamics 
only energy differences can be measured. However, in reacting gases, 
consistent choices must be made, so that the differences in ED values 
between the reactants and the products gives the correct reaction energy 
(at absolute zero). For example, if ED is taken to be zero for rtumic N, 
the corresponding quantity for N* must equal the ionization eneryy of N , 
while that for No must be the negative of the dissociation energy. In 
thermochemistry it is conventional to choose ED = 0 for the elements in 
the phases that are stable at room temperature and pressure. For air - 
(including carbon dioxide) these are gaseous N, ‘ O, . and ar 
and crystalline graphite. The use of graphite as a reference material 
has the disadvantage that it makes the internal energy of gaseous CO, 
negative, and this causes the calculated energies for air near room 
temperature to deviate appreciably from direct proportionality to the 


temperature. To prevent this behavior, ED for gaseous co, (as 


weli as gaseous No P 0, . and Ar) has been set equal to zero in the 
Present work. 
The corresponding values of ED for the atomic species of present 


interest aie given in the first table of the second chapter, toge..ier 


with the dissociation and ionization energies from which they were derived, 
and the corresponding references. An effort was made to obtain the best 


values presently available, since results of thermodynamic calculations 


cone 


are usually more sensitive to inaccuracies in these values than to any 


other inaccuracies. 
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The ideal-gas entropy (at one atmosphere pressure) can be obtained, 
except for a constant term, from the thermodynamic relation 5S = fanr 
(p constent), where H=E+RT and E is given by Eq. (1.3-4). The 


result can be written 
8°/r = (E°- ED)/RT + in Q + : inT+C, (1.3-6) 


where statistical mechanics (Mayer and Mayer, 1940) shows that the 


integration constant is 
a 2a od 
C= 1 + Ink/p® + > in(27M k/2h“) = 2 inM- 2.66496 , (1.3-7) 


with p°= 1 atmosphere and M_ the atomic weight. Other ideal-gas 
thermodynamic properties can be obtained from Eq. (1.3-4) and Eq. (1.3-6) 
using well-known thermodynamic identities. For example, the dimensionless 


Gibbs free energy is given by 


o 


PEED. 50 _ EO Ep s 2 
-l=InQ+t 2 InT + 2 InM - 3.66496 . (1.3-8) 


RT “R RT 


Substitution of this expression in Eq. (1.2-11) gives an equation for the 


equilibrium constant in terms of the partition functions of the reactants: 


3/2 


2 ° ° ° 
N Qh Map Qap E, (AB)~E, (A)-E,(B) 
N f =@\ IM, M 6,0, SF 15 en ak oa (1.3-9) 
AB AB AB 
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When A represents a positive ion, B an electron (with Qn” 2),. and 
AB a neutral atom, Eq. (1.3-9) yields the familiar Saha equation for 
equilibrium ionization. 

In calculating ideal-gas functions by the above equations, however, 
there is a basic difficulty not mentioned in elementary texts. An isolated 
atom or positive ion (in an unlimited volume) actually has an infinite 
number of stable electronic states with energies E; less than its 
fonization energy, and thus the sums in Eq. (1.3-2) and Eq. (1.3-4) diverge. 
However, at temperatures up to several thousand degrees the series are 
semi-convergent in the sense that essentially the same answer is obtained 
by taking the first ten or the first ten thousand terms. The omission 
of the infinite "tail" of the series, whicn contains terms corresponding 
to electrons in highly~excited levels, is justified because in actual 
situations no atom can occupy an unlimited volume, since it is limited 
at least to the volume of the system, and generally, in systems containing 
many atoms, to a much smaller volume, comparable to the mean volume per 
atom (see Section 1.7). 

At somewhat higher temperatures the difficulty is less trivial, 
because there may be no point at which the terms in the sum become 
small, and thus the answer may depend significantly upon where the 
series is brokan off. This is illustrated in Fig. 1.2, which shows the 
Percentage increase in internal energy and free energy as the cutoff point 
is raised from an electron quantum number of n= 4 to n#=8. Various 
workers have suggested simple cutoffs which depend upon the atom density, 


the electron density, or the temperature. However, none of these simple 
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one-parameter relations can be completely correct, because the high 
electronic levels are affected both by close-in neutral particles and, 

to a greater extent, by close-in charged particles, while the temperature 
influences how often these particles come close-in. The true cutoff is 
thus a function of the composition, density, and temperature, so that if 
the sums in Eqs. (1.3-2) and (1.3-4) are sensitive to the cutoff the "ideal- 
gas" functions also depend upon the compositicn and density as wel! as 
temperature, which greatly complicates thermodynamic calculations. 
Moreover, the higher bound levels below the cutoff may be perturbed 
sufficiently to further affect the thermodynamic properties. 

As shown by Fig. 1.2, at successively higher stages of ionization 
an atom'’s properties remain independent of the cutoff to successively 
higher temperatures. in air below standard (sea level) density, as the 
temperature is raised to the point where the cutoff for any particular ion 
becomes {mportant the concentration of that ion becomes so small (due 
to further ionization) that the choice of cutoff has little effect on the 
thermodynamic properties of the equilibrium mixture. Unfortunately, 
this simplification does not hold for the highest air densities of present 
interest (up to 10 times the standard density, as produced by a strong 
shock wave at low altitudes). Even at such densities, however, electronic 
states with ns 4 should be little affected by neighboring atoms, since 
the electronic orbits of such states lie well within the mean interatomic 
or interionic distance. Thus, thermodynamic functions calculated by 
summing up to n= 4 _ represent lower limits to the true values. In 


addition, they usually represent useful approximations to the actual 
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values, since at high densities the correct cutoff is not much greater than 
4 (see Section 1.7), while at low densities the equilibrium properties are 
relatively insensitive to the cutoff, as already pointed out. 

If better, variable~cutoff thermodynamic properties are required they 
may be obtained by adding the contributions from states betwean n= 5§ 
and the cutoff n, - For light atoms and ions, the most important of such 
states are those having only one electron with n> 4 , such states with 
one electron outside a "core" are quite hydrogen-like, with a statistical 


2 times that of the core (with the outer electron re- 


weight equal to 2n 
moved), and an energy just Ry (Z+ 1)2/n2 less than the energy of the 
core, where Ry is the Rydberg (109,737 cm™/ © 13.605 eV), and Z is the 
charge of the atom or ion. The partition-function summation, Eq. (1.3-2), 
over these states can then be factored into a sum over the states of the 


outer electron and a sum over the states of the core. The final result is 


n 
z 
Qh) = Qns4 x) + Qns4 x") exp (- wes an? exp (sgen") (1 237 10) 


n kT 
n=5 


where Q,< Qe) is the partition function for the next higher ion and 

IQ) isthe X- x" ionization energy. Since the cutoff n x depends upon 
temperature, density, and somewhat on composition (see Section 1.7), it is 
more convenient to calculate the correction term on the right hand side of 

Eq. (1.3-10) during the course of an equilibrium ionization calculation than 
to tabulate the "true“ ideal-gas thermodynamic properties of each species 


as functions of two or three variables. 
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Accordingly, ideal gas functions for the atoms and ions of C, N , 
O, and Ar have been calculated for temperatures of present interest 
using a fixed cutoff just above n@= 4. Extensive numerical results are 
presented in Chapter 2. Tables showing the energy~level data used, 
and the fractional population of each electronic state versus temperature 
are also included. The thermodynamic results have been compared with 
results of several previous workers (JANAF Thermochemical Tables, 1960; 
Gilmore, 1955; Woolley, 1957; Kilsk’, Gilmer, and Gilles, 1957; 
Green, Poland, and Margrave, 1960; Yungman et al., 1961; Martinez, 1961; 
Gurvich et al., 1962; McBride et al., 1963; Hilsenrath, Messina, and 
Evans, . 34). Up to 10,000°K, almost all ths results agres to within 0.03 
Per cent. Moreover, the small differences are due primarily to use of 
different values for the gas constant R; if this is corrected for, the 
deviations are generally reduced to less than 0.01 per cent. An exception 
is the tabla for “I! by Hilsenrath, Messina, and Evans (1964), which gives 
values deviating from other work by up to 2 per cent (their values for the 
other air atoms and ions show no such deviation). 

Only a fev «ples extending above 10,000°K have been published 
previously, and these show increasing variations in the results for the 
neutral atoms, depending upon the electronic cutoff used by the different 
investigators.( Variations fcr the ions appear only at considerably higher 
temperatures, as shown by Fig. 1.2.) The earlier work by the writer (Gilmore, 
1855) up to 24,000°K used a cutoff around n= § but omitted some states. 
Up to 12,000°K those results agree very well with the present results, but 


differ by up to 6 per cent for the neutral atoms at 24,000°K. The values of 


Woolley (1957) to 4,000,000°K and the extensive Russian tables (Yungman et al., 
1961; Gurvich et al., 1962) to 20,000°K also show deviations above 10,000°K. 


Woolley (1957) used a somewhat differant type of cutoff that excludes many 
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states included in the present work (planning to include them in subsequent 
equilibrium calculations). Hence, his values for the enthalpy, ontropy, 
and negative free energy are generally lower than the present values; for 
the neutral atoms the difference is a few per cent or less at 20,000°K. 

On the other hand, the Russian values aro generally higher because they 
used cutoffs at n= 11-13 , based on restricting the electrons to the 
mean volume per atom at densities corresponding to a pressure of one 
atmosphere. At 20,000°K their enthalpy values for the neutral atoms 

are almost a factor of two higher than the present values, although their 
entropies and free energies are only a few per cent higher. Whether such 
differences are significant in equilibrium thermodynamic calculations 
depends upon the density under consideration and also upon whether separate 
allowance for the highly excited states is made in the equilibrium computa- 


tion, as already discussed (see also Section 1.7), 


1.4 Ideal-gas properties for diatomic gases 

The relations Eqs. (1.374) and (1.3-6) for the energy and entropy of 
a monatomic gas apply equally well to a molecular gas, provided that the 
summations are extended over all the molecular energy levels, which can 
differ ‘not only in electronic energy but also in rotational and vibrational 
energy. Because of these additional degrees of freedom, the individual 
levels of a molecule are so numerous that it is rather impractical to 
tabulate all of them. Fortunately, however, for each degree of electronic 
excitation the rotational and vibrational levels are usually quite regular 


and can be represented by simple formulas, only the coefficients of which 
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need to be tabulated. Moreover, using these formulas, the terms of the 
series in Eqs. (1.3-2), (1.3-4) and (1.3-6) can often be summed 
algebraically with reasonable accuracy, at least in some temperature ranges. 
The energy levels ct molecules can be grouped into electronic states, 
each of which is characterize (in part) by the total spin 8 of the 
elestrons. ‘vx" as for vioms, the multiplicity 28% 1 , is written as a 
superscript on the left of the state symbol. Unlike the atoms, however, 
diatomic molecules have an internuclear axis, and the second important 
quantum number is not the total orbital angular momentum, L , of the 
electrons, but the component , A , of this momentum in the axial direction. 
In analogy with the Roman-letter designations of atomic terms, diatomic 
states are designated by the Greek letters ©, ,4, 8, .+., corresponding 
to A@0,1,2,3, «+++ Because in diatomic molecules terms of a given 
L are separated into states of different A , the statistical weights of 
these states are generally less than those of the corresponding atomic 
terms; specifically, the weights are 28 +1 for <= states and 2(28 + 1) 
for 1 , A ,#, «-, states. States with S>0 and A>O_ split 
into 28 + 1 substates, designated by writing the vector sum of A 
and the axial component of S as asubscript. However, this splitting 
can often be ignored, since it is usually small and often associ ated only 
with higher excited states. Other subscripts and superscripts are some- 
times added to the state symbol to indicate the electronic symmetry, but 


this does not affect the present considerations. 


22 


Bince a diatomic molecule may have more than one state with © 
the same 8 and A values, and also since there may be a considerable 
period between the experimental discovery of & state and the dutermi nation 
of its type, the states are also labelled somewhat arbitrarily by Roman 
letters, with X designating the ground at"te, and A, B, GC, ... the 
successively higher (or earlier discovered) states of the same multiplicity 
as the ground state, while a, b,c, ... designate states of a different 
multiplicity. (Exceptions are the N. and C, mciecules, where 
the capital and lower case letters are reversed, due to early misidentifi ca~ 
tion of the ground state.) 

The electronic energy of a diatomic molecule depends not only upon 
the electronic state but aleo upon the interatomic (or, more accurately, 
internuclear) distances. Curves showing this vari ation for N, and ND 
are presented in Fig. 1.3. At large internuclear distances the energy must 
be that of two individual atoms. At distances approaching zero the energy 
must approach infinity because of the strong Coulomb (electrostatic) 
Tepulsion between the two nuclei. At intermediate distances the curve 
must have a minimum if the molecular sta‘e is to be stable. The two nuclei 
will tend to approach this minimum-energy distance, but may vibrate about 
this point. Since the vibrational motion of the nuclei is much slower than 
the orbital motions of the lighter electrons, the electrons will keep adjusted 
to the instantaneous position of the nuclei. Thus, the electronic energy 
will follow the same curve regardless of the amplitude of vibration, and 


such a curve forms an effective potential for molecular vibration. 
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The lower portion of an attractive potential curve (see Fig. 1.3) 
can often be approximated by a parabola. Quantum theory (Herzberg, 1950) 
shows that the vibrational energy in a parabolic potential equals w(v + 1/2), 
where is @ constant which depends inversely upon the width of the 
parabole, v is the vibrational quantum number, which can take on the 
values 0, 1, 2, ..., and the sero of energy is measured from the bottom 
of the parabola. Actual potential curves, however, diverge more and more 
from a parabola at higher energies; therefore the vibrational energy levels, 
instead of being cvenly spaced as in the above formula, fall closer and 
closer together as the dissociation limit is approached. It is conventional 
to fit these levels by the formula 


Glv) = wv 1/2) -wxlv+ 1/2)? + ay (vt 1/2)74..., 
(1.4-1) 
where wee WX 7 Ww ava?! etc., ate constants determined spectroscopically, 
and the subscripted e indicates that this is an expansion about the 
equilibrium point (potential curve minimum). The short numbered lines on 
the potential curves of Fig. 1.3 indicate the observed vibrational energy levels 
of these molecules. 

Besides vibrating, a molecule can also rotate. The simplest approxi- 
mation is that of a “rigid rotator, " whose quantum-mechanical energy levels 
(Herzberg, 1950)are given by BJ(J+1), where B iz a constant inversely 
proportional to the square Of the internuclear distance, and the rotational 
quantum number J equals 0, 1, 2, ... Since actual molecules are 


not rigid, but stretch as they rotate so that the effective value of 8B 
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decreases with increasing J , a negative correction tarm proportional 
to ry > 1)? is conventionally added to the energy expression. Moreover, 
if the molecule is also vibrating, this affects the mean internuclear 
distance and gives a B value which depends somewhat on the vibrational 
quantum number. 

By adding the contributions of electronic, vibrational, and rota- 
tional energy one obtains the following general] formula for the molecular 


energy levela (term values): 


on 


A (vi =o + a, (vtI/2) ~ wx, vei/2)? + wy (ve1/2) ® 
(4-2) 
t..0% BIH) - D704? + HPO + ..., 


where 


Bo = BL - a (v+l/2) + ¥e(vti/2)? $e) 


= / 
D, D, + B (vt 1/2) + ..., 


Hy ™ H, + 6, (v+ 1/2) + «s- ’ 


and “> : is the electronic energy at the equilibrium distance, taking the 
zero of ensrgy to be the lowest level (v= 0, J = 0) of the ground state. (The ~ 
asterisk is added to distinguish this quantity from the conventional ie a4 
the energy above the potential minimum of the ground state.) In Eq. (1.4-2) 


some generally-negligible terms due to the interaction of the electronic and 
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rotational angular momenta have bee. omitted. 

The statistical weight of each level is just the electronic statistical 
weight already discussed, multiplied by 2+ 1 for most molecules. 
However, if the ~olecule is homonuclear (like Ny : Ny , otc.) the 2+ 1 
is replaced by a rapidly-oseillating function of J , which for present 
purposes can be approximated by its mean value, (2) + 1)/2 . 

Siace each electronic state has only a finite number of bound 
rotational-vibrational levels, : the summations required for thermodynamic 
calculations can be carried out without eny further convergence difficulties. 
However, at temperatures cf a few thousand degrees, many hundreds of 
terms make significant contributions to the sums, so that hand computations 
become lengthy. For this reason, Mayer and Mayer (1940) have worked 
out approximate algebraic formulas for these sums, based on replacing the 
summation over rotational leveis by an integration, and neglecting or 
approximating the higher correction terms in the energy~level formulas. 
However, some of these approximations become poor at very high temperatures. 
Moreover, with modern high-speed computers it is virtually as easy to perform 
the summation directly. 

At temperatures above 5000°K or §000°K, vibration-rotation energy 
levels of the lower electronic states near their dissociation limits can 
make significant centributions to the thermodynamic properties. Unfortunately, 
spectroscopic measurements on such levels have n>t been made for most 
electronic states of interest. For the few states for which they are 
avaliabic, it is found that a large number of terms must be included in 
Eq. (1.4-2) in order to fit the data. This indicates that the usual measure- 


ments on only the lowe- levels of an electronic state cannot be safely 


Neglecting a few highly~excited states which dissociate to positive 
plus negative atomic ions. 26 
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extrapolated to get the very high levels. However, these levels are 


determined by the high portions of the corresponding potential curve. 
It is often possible to determine this curve with some confidence, by 
use of Rydberg~Klein and valence-bond calculations (Vanderslice, 
Mason, and Lippincott, 1959; Gilmore, 1965). The energy levels 

| corresponding to this curve can then be calculated quantum- 
mechanically. 

For present purposes, however, it 1s possible to bypass this last 
step, and calculate the contributions of the high levels to the thermodynamic 
properties directly from the potential curve. At temperatures where these 
levels contribute significantly, their spacing is small compared to kT, so 
that the formulas of classical statistical mechanics can be used. These 
formulas show that the partition function for structureless particles of mass 
m with total energies between Ei and EL + inavolume dV where 


the potential is U , is given by(Mayer and Mayer , 1940): 


where erf is the error function. 


diatomic energy levels between E, 


The partition function contribution from 


and the dissociation energy D can 


be obtained by multiplying the above expression by the electronic statistical 


weight Ga letting m_ be the rediced mass, M,M,AM,t+M)Q 


27 


(1.4-3) 


E 
2 
3/2 
em) | (E~W) 1/2 exp (-E/kT) dE 
h 

Ey 

(E,-U)AT 
2ntamen)?/ey eo WkT [¢ erf x2/2 . y 142 | ' 

h (E,-WAT 


lien et 


wage thee 


Ree 


Sites OP oe a nile —— ~<a if 


one ee 


atrees. 


ee ce ere Raa A ot EGE tT eS Mp ores. 


where M, and M, are the atomic weights of the two atoms, replacing 


dV by anc*dr where r is the internuclear distance, and integrating 


over r. The result is 


= r 
Qe, to D) = 5 fn()r dr - f, (r)r ar], (1.4-4) 
To % 


where 


M2 -vtdkT _ V2 pt 
kT 


£, (r) is the game function except that D is replaced by E, + T> is 
the point on the inner branch of the potential curve where U(r)= D , and 
TY) and rT. are the two points where U(r) = Ey - Asimilar but slightly 


longer expression can be derived for the summation appearing in Eq. (1.3-4). 


The above approach omits the quasistable rotational levels above 
the dissociation energy (Herzberg, 1950), which seme investigators include. 
However, such levels are more conveniently treated as two separate atoms 
subject to an interatomic attraction (see Section 1.6). Moreover, there 
is little point in including bound levels near or above the dissociation 
energy unless the unbound or “repulsive" states or levels of similar energy 
and internuclear distance are also included, since the two types make 
comparable contributions to the equilibrium thermodynamic properties. 


Inclusivn of the latter type, however, requires a departure from the ideal- 


gaa approximation. 
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To check the accuracy of the above relations, partition function 
and energy calculations for several of the states of N» , NO , and O, 
were made using the direct summation method up to 1/3 , 1/2 , and 2/3 
of the dissociation energy, respectively, and the classical integrals 
beyond these energies. The results using the three different crossover 
points were virtually identical up to several thousand degrees, and agreed 
within 0.5 per cent for the entire range from 1000 to 40,000°K. 

The total partition function end energy of diatomic moiecules may 
be obtained by adding such contributions from all the electronic states. 
Of course, just as for atoms, isolated molecules have an infinite number 
of high-excited electronic states, while at finite densities these are 
“cut off" in some complicated fashion by electron-ion interactions. 
However, for most molecules such states lie considerably above the 
dissociation energy, so the molecules tend to dissociate before the 
choice of electronic cutoff makes much difference. In fact, unless the 
repulsive as well as attractive states near the dissociation energy are 
treated carefully, there is little point in including states above the 
dissociation energy (unless their fractional populations are desired for 
other purposes, such as radiation calculations). 

, In the present work, thermodynamic calculations were made for the 
diatomic molecules Ny, Nj, NO, NOY 0}, Oy, Of, and CO, which 
make a significant contribution either to the thermodynamic properties 
or to the charged-particle concentration of equiiibrium air. These 


calculations included all the known and predicted bound states up to the 


29 


aed eneeellciatiinkoninnnemeer pemennntains ieeeeteeraniandioin conte a 


lowest dissociation limit, and in some cases one or two states above | 
this limit. The method used was to sum over the enercy level up to half 
the dissociation energy of each state, where the levels were calculated 


using the spectroscopic data given in Table 86 of the supplementary 


volume. (The higher rotational constants D, and H,, for which 
data are sparse, ware calculated from the vibrational constants using 
gormulas given in Herzberg (1950).) Levels above half the dissociation 


energy were included by means of the classical integral, Eq. (1.4-3), 


using the potential curves of Figs. 1-3 to 1-5, and of Krupenie and 
Weissman (1965) for CO. 

The results are presented in Chapter 2. Up to 6000°K they 
generally agree very well with those of previous workers (Gilmore, 
1955; Beckett and Haar (1957); JANAF Thermochemical Panel 
(1960); Yungman et al., 1961; Gurvich et al., 1962; Bristow and 


McChesney, 1965) except for small differences due to use of more recent 


values for the gas constant R and the second-radiation constant Cy = he/k . 

Above 6000°K, however, the present values for most molecules begin to 

diverge from the older values, usually in the positive direction due to 

inclusion of more electronic states in the present calculation than in any 

previous calculation, except that of Bristow and McChesney (1965). At 

still higher temperatures, around 15,000 or 20,000°K, the present values 

cross over and fall below the few previous values available, because 

previous investigators effectively included some levels above the dis- 


sociation limit of each state. 
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1.5 Ideal-gas properties for polyatomic gases 

The only polyatomic molecules which contribute as much as 0.01 
per cent to the equilibrium thermodynamic properties of dry air are co, 
and NO, . In addition, the negative ion NO) can affect the electron 
concentration and hence the radio-wave absorption of high density air 
at a few thousand degrees, so it is desirable to include it in the 
equilibrium calculation, 

The ideal-gas properties of polyatomic molecules can be computed 
by summing over their rotational, vibrational, and electronic energy 
levels in the same way as already described for diatomic molecules. 

Of course, the additional rotational and vibrational degree of freedom 
of the larger molecules produce a more complex set Of energy levels. 
However, since the concentrations of polyatomic molecules generally 
become small above 5000 or 6000°K, the calculations can be restricted 
to lower temperatures, where excited electronic states and high rotational 
and vibrational levels can usually be neglected. This permits relatively 
simple approximations to the ideal-gas thermodynamic properties 
(Mayer and Mayer, 194C). 

Results of such calculations are available in the JANAF Tables 
(1960) for the polyatomic molecules of present interest except NO, . 


Thermodynamic values for the latter are tabulated in a recent report by 
Clifton (1966). 
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Since the existing ideal-gas tables for polyatomic molecules are 
fairly adequate for high-temperature air calculations, no further 
equations or tables will be presented here. 


1.6 Effects of interparticle forces on the thermodynamic properties of air 

In the ideal-gas approximation the interactions between the 
molecules, atoms, ions, and electrons in high-temperature air are 
neglected, except when two or more particles are bound together and 
can be treated as a single particle. This approximation is reasonable 
at low air densities, but at sufficiently high densities the mean distance 
between free particles becomes so small that such interactions can no 
longer be neglected. For neutral molecules and atoms, interaction 
forces are very small except at distances less than about twice the 
intermolecular distance in the liquid or solid phase. The highest air 
density associated with nuclear fireballs or missile flow-fields is that 
produced by a strong shock wave at sea level, which is roughly 0.01 g/cm? 
(about 10 times the ambient density). Since this density is about 1 per cent 
of that of liquid air, one may expect that neglect of intermolecular forces 
will produce errors of the order of one per cent in the thermodynamic 
properties. 

At temperatures s0 high that the air is largely fonized, considerably 
greater errors can be made by neglecting the interactions between the 
fons and electrons. This large effect is due to the long range of the 
Coulomb forces between cherged particles, which decrease with distance 
Uke 1/r? , in contrast to the forces between neutral molecules, which 


decrease like ve! for large values of r. As an example, consider 
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air at 0.01 g/cm” and 70,000°K, where the major equilibrium species 
are n* ; o* , and free electrons. The ideal-gas pressure under 

these conditions is about 10,000 atm, while the total Coulomb force 
between neighboring electrons and positive ions (at their mean distance 
On opposite sides of a unit surface) is about 1000 atm. Accordingly, one 
may axpect errors of the order of ten per cent in the ideal-gas approxi- 
mation under these conditions. This makes it quite desirable to include 
corrections for charged-particle interactions in the thermodynamic 
calculations for air, as well as somewhat desirable to include neutral- 
particle interactions. 

The standard method for the thermodynamic treatment of moderately 
dense gases is to write a virial expansion, where the first term is the 
ideal-gas contribution and subsequent terms give the contributions from 
two-particle, three-particle, etc., interactions. The virial equation 
for the pressure of a gas mixture may be written (Mayer and Mayer, 1940; 
Hirschfelder, Curtis, and Bird, 1954) 


- pe NRT p> NN,B,,(1) RT + LLL RNA,C,,,(1) RT+..., (166-1) 


1 


where By , (7) ’ C55, (1) 1 eee are called the second, third, etc., virial 
coefficients. These coefficients can be evaluated in terms of the inter- 
particle forces. If we treat the air particles as spherical so that they 
have a two-particle interaction energy Uy which depends only on 


their distance r , the second virial coefficient is given by (Mayer and 
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Mayer, 1940; Hirschfelder, Curtis, and Bird, 1954) 
B, , (7) = ma fh - ex9(-U,,(01/47)] adr, 
° 


where a 4s Avogadro's number. The expressions for the higher virial 
coefficients are much more complicated, but fortunately the corresponding 
terms in the virial expansion can generally be neglected for air densities 
of present interest. 

A virial expansion can also be written for the internal energy, similar 
to Eq. (2.6-1) for the pressure, except that temperature derivatives of the 
virial coefficients, TAB, ,/aT . etc., appear instead of the coefficients 
themselves. Tabulations of virial coefficients for various molecules 
(Hirschfelder, Curtis, and Bird, 1954; Woolley, 1962) show that, at the 
high temperatures of present interest, such derivatives are typically an 
order of magnitude smaller than the coefficients themselves. This is to 
be expected, since the longer-range intermolecular interactions {except for 
Coulomb interactions) are fairly small compared to thermal : -rergies at such 
temperatures, so that molecules behave roughly like rigid spheres, with 
virial coefficients approximately independent of temperature. Accordingly, 
for air densities and temperatures of present interest, it is reasonable to 
neglect virial corrections to the internal energy (except for charged-particle 
interactions; see Section 1.7). 


The Helmholz free energy is given by 


A= Aeal * vaT OX N,N, B, (7) ate 3p (1.6-3) 


mixture 
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where V_ is the volume of the system. The virial corrections to the 
entropy and the Gibbs free energy may be obtained readily from Eqs. 
(1.6-1) and (1.6-3) with the help of the thermodynamic identities 
S=(E-A)/T and F=«A+ pV. 

In Section 1.2 the ideal-gas expression for the chemical-equilibrium 
constant, in terms of partial pressures, was derived by minimizing the 
Gibbs free energy. For nonideal gas mixtures, however, a similar deriva-~ 
tion is not convenient because the total pressure is no longer the sum of 
the individual partial pressures. The equilibrium constant may be obtained 
instead, by minimizing the Helmholz free energy at constant temperature 
and volume. Equation (1.6-3), after substiti:tion for the ideal free energy 


from Eq. (1.2-4), can be written 
A/VRT ps N, [rover - 1+ In 8,20] p> RNB, (1) . (1.6-4) 


The minimization condition that the differential vanish yields, after seme 


cancellation, 


0-5 [Fever + InGi,RT) + > Nya] oN, (1.6-5) 


For variations due to a single chemical reaction, XY“ X+Y , one has 


oNyy* - éNy = - oN, » while all other 5N, vanish, so that Eq. (1.6-5) 


yields 
N Fo_-Fo- Fe a . 
a = RT exp |- Sia Ta - 2 5, (Ber, s-®"®) ; (1.6-6) 
Y 


The ideal-gas part of this equation agrees with Eq. (1.2-11). 
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In order to apply the relations derived above to thermodynamic 
calculations, numatical vaiues for the second virial coefficients, 

B,;(T), are needed. For most common molecules that are stable at 
room temperature, values of the virial coefficients for like molecules 

(i = j) and a few for unlike molecules (i # j) have been measured over 
the easily-accessible temperature range, and extrapolated to higher 
temperatures by fitting an intermolecular potential according to Eq. 
(1.6-2) (Hirschfelder, Curtis, and Bird, 1954). Woolley (1962) has 
calculated and tabulated such values for air molecules up to 15,000°K. 
However, at high temperatures such extrapolations often give too high 
values. Better virial coefficients for N, and the rare gases have been 
calculated by Amdur and Mason (1958), using potentials derived from 
molecular beam scattering. In the calculation range of 1000 to 15,000°K 
their values for No range from 0 to 19 per cant lower than those of 
Woolley, while their values for Ar are 10 to 23 per cent lower. 

In high temperature air calculations, virial coefficients are also 
needed for species such as atomic N and O , for which no measure- 
ments are available. Woolley (1962) has also made estimates of these 
coefficients by deducing the interatomic porentials from the corresponding 
intermolecular potentials. For interactions noi involving chemical bending, 
suchas N- No or O-Ar , his results are not unreasonable. However, 
better values for the N ~- No virial coefficient may be obtained by using 
the potential calculated by Meador (1960) using valence-bond theory. The 
results range from 25 to 40 per cent lower than Woolley's results, over the 


temperature range from 8000 to 15,009°K. It may be added that Meadorc's 
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potential, and the derived virial ccefficient, are probably about 4 per 
cent too low, judging by a comparison of his No - No potential with 
the experimenta! results of Amdur, Mason, and Jordan (1957), 

Nonbonding potentials, at temperatures of several thousand degrees 
or more where the small van der Waals attraction is negligible, can 


usually be fit by an exponential repulsion: 


U(r) = ae , (1.6-7) 


where a and c are constants which depend upon the two molecules 


or atoms. ‘She integral in Eq. (1.6-2) can then be approximated by (Amdur 


and Mason, 1958) 


an 
B(T) = << ing 4810] , kTK<<a. (1.6-8) 
Cc 


The recommended values for the constants, obtained by increasing 
Meador's (1960) values by 4 per cent, are 

N,-N,: a=8 s= 2.78 x 10° cm! 

27 Ny: 8 33 eV, 2 2.78 x cm : 


N -N,: a= 363 eV, c=2.85x 108° cm? . (1.6-9) 


Much less is known about the high-temperature interactions of 
other air molecules and atoms, although some theoretical calculations 
involving O, . NO , and O are available (Meador, 1960). 
However, these molecules and atoms, and even many of the minor air 
species like CO and C _ , are approximately the same "size" as 


No or N_ , and should have roughly the same high-temperature 
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virial coefficients. For air above 8000°K, it is reasonable, then, to 
use Eqs. (1.6-8) and (1.6-9a) for all neutral molecule-molecule inter- 
actions, and Eqs. (1.6-8) and (1.6-9b) for all neutral atom-molecule 
interactions. 

Virial coefficients for interactions between neutral and charged 
Particles, on the other hand, are generally much smaller (Woolley, 
1960) because the polarization attraction counteracts the core repulsion. 
For present purposes it seems adequate to set them equal to zero. 
Coulomb interactions between charged particles will be considered 
later, in Section 1.7. 

Interactions involving chemical bonding, such as the important . 
N-N ,N- OO, and O-O interactions, cannot be adequately | 
treated by considering only a single interaction potential. Instead, as 


shown in Figs. 1-3, 1-4, and 1-5, the two atoms may approach each 


| 
«her on any of several potential curves, depending on the relative | 
3 orientation of their electrons. The correct virial coefficient to use is | | 
{ a weighted average of the values computed using the various curves, 
where the proper weights are the statistical weights of the molecular 

states (see Section 1.4), and contributions of those states or levels 
: already included in the molecular partition function should be omitted 
here. 

The present calculations of the partition function and thermo- 
dynamic properties cf diatomic molecules, described in Section 1.4, 
include for all the lower electronic states every rotational-vibrational 
leve! below the corresponding dissociation energy. For attractive | 


potentials, to be consistent, the virial integral of Eq. (1.6-2) must then 
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be replaced by one which omits contributions from bound states. The 
result, after taking the weighted average over the different electronic 


states n , is 


wo = /kT is 
z nS ee eee _ay 1/2 9-E/kT 
B(T) go - nl * dr j1-fe Lae (E-U ve dE 
fa U,70 U< 0 
(1 26-15) 


where for simplicity the subscripts ij designating the interacting 
atoms, and the variation of UL with r , are not explicitly indicated. 
The last integral in Eq. (1.6-10) can also be expressed in teri.t of the 
error function, as in Eq. (1.4-3). 

Sample calculations for nitrogen and oxygen aioms at high tempera- 
tures were carried out using Eq. (}.6-10). They gave second vir‘al 
coefficients almost an order of magnitude smaller than those involviag 
molecules, because negative contributions from nonbound levels in the 


attractive potentials largely cancelled the positive contributions from the 


repulsive cores. Accordingly, for present purposes the atom-atom virial 
corrections can be neglected. 

It might be mentioned that when high-temperature diatomic partition 
functions are calculated by integrating cver all portions of the potential 
curves and all energies, as done by Beckett and Haar (1957) and Bristow : 
and McChesney (1965), instead of cutting tham off at the dissociation 
energy as done in Section 1.4, the corresponding atom-atom virial correction 
is already included implicitly in the molecular thermodynamics. The only 


remaining contributions to the coefficient are those from any electronic 
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states not included in the diatomic integration. In principle, the 
diatomic Partition functions could include alj electronic states, thus 
making the atom-atom virial coefficients vanish identically, but in 
practice the purely-repulsive states are usually omitted from the 
molecular treatment. The corresponding second virial coefficients, 
obtained by summing Eq. (1.6-19) only over the repulsive states, are 
somewhat larger than those obtained by using the complete sum, but 
still only about half as large as the coefficients estimated by Woolley 
(1960) by scaling down the molecule-molecule potentials. 


1.7 Effects of Coulomb forces on the thermodynamic properties of 
ionized air 


The Coulomb forces between the ions and electrons present in 
high-temperature air are sufficiently different from the intermolecular 
forces already considered as to require a separate treatment. When 
U,,) = constant/r , the integral in Eq. (1.6-2) is found to diverge at 
r=0Q when the constant is negative (i.e., for charged particles of 
opposite sign), while for large r it behaves like +t f rdr and diverges 
for both positive and negative Coulomb potentials. The divergence at 
r=0 is due to inclusion of bound states; it may be removed by including 
only the states with energies above that of the separated particles, by use 


of Eq. (1.6-10). The divergence at r=” , however, is more funda- 


mental. A uniform gas of charged particles which ‘s not electrically neutral 


(i.e., has more positive than negative charges, or vice versa) can be 
shown to have a Coulomb energy per unit mass which depends upon the 


size and shape of the gas volume considered, and becomes infinite as 
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the volume goes to infinity (at fixed gas density). It is not curprising, 
then, that the virial coefficients for a Coulomb potential diverge. In an 
ionized gas with no net charge, the infinite positive and negative virial 
terms from the repulsive and attractive Coulomb forces, respectively, 
must somehow cancel to first order, leaving only a finite remainder. 

If the charged~particle density is not too high, this problem can 
be treated by the Debye-Huckel theory, which determines the mean 
distribution of electrons and ions around any given electron or ion, using 
a linearized self-consistent-field approximation. Since the derivation 
and results are available in several texts (Fowler, 1936; Fowler and 
Guggenheim, 1956; Cambel, Duclos, and Anderson, 1963) they will not 
be reproduced here. Qualitatively, the effect of Coulomb interactions is 
to decrease the effective ionization potential of each species, and thus 
increase the degree of ionization. The interactions also decrease the 
pressure and energy of the ionized gas mixture directly, in contrast to 
the indirect increase due to the increased ionization, so that the net 
thermodynamic corrections are positive at some temperatures and densities 
and negative at others. 

For present purposes, however, the conventional Debye-Huckel 
treatment needs to be supplemented by two additiona! considerations. 
The first concerns the dielectric constant. The original analysis of Debye 
and Huckel was applied to ions in solution, and since the force between 
adjacent ions is affected by the polarizability of the intervening liquid, 
the dielectric constant of the latter enters into the basic equations. 

A similar effect may be expected in an ionized gas whenever molecules, 


atoms, or polarizable ions are situated between neighboring ions, but in 
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gases the continuum fluid approximation and the use of the static 
dielectric constant will be less accurate. Fortunately, at gas densities 
and temperatures of present interest, the dielectric constant will not 
deviate more than a couple per cent from that of empty space, so that 
it can be approximated by the latter value, with an error less than the 
other errors inherent in the Debye-Huckel approach. 

The second consideration involves the electronic-state "cutoff" 
introduced in Section 1.3. Conventional Debye-Huckel theory does not 
consider electronically-excited states. It does, however, give fer the 


mean potential about each ion of charge Z : 


U(r) = (Ze/r} exp (-r/d) , (1. 7-1) 


where d is the Debye length, given by 
d? = ktT/ (sne*5 nz?) : ; (1. 7-2) 


and the summation is taxen over all charged particles, including electrons. 
Some workers have assumed that the higher excited states of each atom or 
ion are precisely the bound states of an electron in the Debye-Huckel 
potential, Eq. (1.7-1). This potential has only a finite number of bound 
states, so that the partition function converges. However, these states 
include a number having classical electron orbits with radii approaching 
the Debye length. Since the Debye-Huckel treatment is valid only when 
the Debye length is longer than the mean distance between ions (Cambel, 
Duclos, and Anderson, 1963), such states have a bound electron which ts 


generally closer to other ions than to the ion to which it is presumed to be 
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bound. This is not reasonable. Moreover, from the viewpoint of 
classical statistical mechanics, in which an electron is characterized 
by a position and a momentum or energy, it is clear that such a procedure 
will give a partition function which counts porsions of phase space more 
than once. 

A better procedure is suggested by the classical approach: The 
entire volume of the gas can be divided into approximately spherical 
volume elements surrounding each ion, and the bound states within each 
element calculated. The volume elements may be taken of equal size, 
but when ions of different charges are present it is preferable to make 
the size of each element proportional to the charge of the ion it contains. 
Since the radii of the volume elements are smaller than the Debye-Huckel 
length d , it is reasonable to omit the exponential factor in Eq. (1.7-), 
leaving just the ordinary Coulomb field. The lower bound states of electrons 
in this field will be just the Rydberg states already described in Section 1.3, 
while the higher ones wili deviate due to the finite volume. However, for 
gas densities of present interest the latter states are high enough to be 
approximated classically. Classically, the bound electron is unaffected 
by the finite volume until its energy is sufficient to permit it to reach the 
boundary, at which point it ceases to be beund. This principle gives 


a cutoff quantum number for substitution in Eq. (}.3-10): 


2 2 = 
A= Ry (Z+ Nro, 7e ’ (1. 7-3) 


where T+} is the radius of the volume element about an ion of charge ! 
Z+1 (Z+ 1 appears instead of Z because a highly-excited atom or ion 


of charge Z is treated like an ion of charge Z+ 1 plus an outer electron.) 
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1.8 Equilibrium calculations and results for air 

Based on the theoretical relations presented in the precediag 
sections, a FORTRAN code has been written and equilibrium calculations 
carried out for air between 10,000 and 19,000,000°K, including the 


Debye-Huckel, variable electronic cutoff, and second virlal corrections. 


The iteration process selected, based on the assumption that tre molecular 


concentrations are smaller than the atomic ones, failed to converge for 
teinperatures below 10,000°K at high densities. No attempt was made to 
modify it, due to lack of time and the adequacy of the earlier results of 
Hilsenrath and Klein (1963, 1965) in this lower temperature region. 

The basic composition of normal] air, used as input, was about the 
game as that used previously (Gilmore, 1955, 1959; Hilsenrath and Klein, 
1963, 1965), except that the CO, concentration was decreased slightly 
in accordance with recent measurements (Kelley and La Chapelle, 1966), 
and the small amount of neon was lumped in with the argon. This gives 


the initial composition shown below: 


Molecule sss SC Mle per cent 


No 78.084 
O, 20.946 
Ar 0.93€ 
CO, 0.032 
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The density range covered was from 1079 


to 10 times standard density, 
where the standard density is 1.2923 x 1073 g/om", corresponding to the 
ideal gas mixture at 273 .15°K and 1 atm pressure, At the lowest density 
considered the assumption of thermodynamic equilibrium might seem unrealistic 
in most practical situations, However, for energy densities corresponding to 
temperatures above about 10,000°K many free electrons are usually present, 
and characteristic times for these electrons to equilibrate with each other 
and with the ion excitation and ionization are only a few milliseconds, Of 
course, equilibration of ion velocities will take a longer time (though often 
still short compared to radiative or hydrodynamic cooling times of large 
masses of air), but the ion kinetic energies make only a relatively smali 
contribution to the total energy and pressure at the higher temperatures. 
Radiation will also depopulate nonmetastable states below their equilibrium 
concentrations, but even in equilibrium their populations are small, since at 
low densities ions tend to ionize further rather than become excited, 

The calculated equilibrium compositions are presented graphically in 
Figs. 2.1 to 2.8. The concentrations of most species agree within a percent 
with the values of Hilsenrath and Klein (1963, 1965) up to 15,000°K, after 
correcting for the 3 percent smaller initial CO, concentration used in the 
present work. There are, however, differences of severai percent in the 
molecular concentrations, due to use of improved diatomic thermodynamic 
properties and virial coefficients. Moreover, the O concentration is 1p 
to a factor of 2 lower, and the OF concentration a couple orders of magnitude 
higher due to use of revised electron affinities. However, at all temperatures 


and densities considered the negative ion densities are more than an order 


+ erect = 


of magnitude beiow the electron densities, 
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The calculated equilibrium thermodynamic properties are presented in 
Tables 102 to 111. These values agree within 1 percent with those of 
Hilsenrath and Klein up to 15,000°K. They also agree within 2 percent with 
the earlier values of the writer (Gilmcre, 1954, 1959) up to 24,000°K, and 
of Hilsenrath, Green, and Beckett (1957) up to 5,000,000°K, except for 
differences of several percent at the higher densities due to omission of 


virial and Debye-Htickel corrections from the earlier work. 
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FIG. 1-1 ENERGY LEVELS OF THE NITROGEN ATOM BELOW 30 eV 
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Chapter 2. THE EQUILIBRIUM THERMODYNAMIC PROPERTIES 
OF HIGH-TEMPERATURE AIR: TABLES AND GRAPHS 


This Chapter contains the following material: 


1, Energy Zero for Air Species: Table l. 
2. Ideal-Gas Properties for Monatomic Air Species: Tables 2-45, 
3. Energy Levels and Fractional Electronic Populations for Monatomic | 
Air Species; Tables 46-85, 
4. Spectroscopic Constants for Diatomic Air Species: Table 86. | 
| 
5. Ideal-Gas Properties for Diatomic Air Species: Tables 87-94. i 
6. Fractional Electronic Populations for Diatomic Air Species: 
Tables 95-101. 
7. Dimensiopless Pyessure, Pressure and Log of Pressure of Equilibrium | 
Air, 10% to 10’ °K 1079 to 10 Fold SeaLevel AirDensity: Tables 102-104, 
8, Dimensionless Internal Energy, Internal Energy and Internal Energy 
Density of Equilibrium Air, Same Temperature and Density Ranye 
As Above: Tables 105=+107, 
9, Enthalpy of Equilibrium Air, Same Temperature and Density Range 
As Above: Table 108. a 
10, Effective Gamma of Equilibrium Air, Same Temperature and Density 
Range As Above; Table 109, 
ll, Dimensionless Entropy of Equilibrium Air, Same Temperature and 
{ 


Density Range as Above: Table 110. 


12. Equilibrium Composition of Air, 104 to 107 °K, 107® fold to Normal 
Sea Level Air Density: Figs. 2-1 to 2-8, 
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Species 


co, 


E2(cm™}) 
0 

3,065 
-29,100 


0 
125,667 
7,506 
82,253 
~3,470 
0 
97,295 
23,351 


92,315 
183,129 
379,788 
766,002 
128,179 

4,448,629 
8,400,690 


Energy of Air Particles at 0°K 


Species 


ES (om™ 4) 
39,359 
156,573 
395,324 
777,950 
1,402,801 
2,192,317 
6,645,117 
12,025,206 


8,705 
20,630 
130, 467 
413,711 
856,518 
1,480,915 
2,399,617 
3,513,617 
9,476,617 
16,505,010 
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Species 


E8(cm™ 4) 

0 

127,110 
349,958 
679,924 
1,162,300 
1,767,400 
2,504,000 
3,504,000 
4,661,800 
8,069,100 
11,930,000 
16,277,000 
21,263,600 
26,797,400 
32,892,900 
39,787,100 
47,196,400 
80,386,400 
116,086,300 
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FIG, 2-1 EQUILIBRIUM COMPOSITION OF AIR AT p/p, = 10-6 
SOLID CURVES; USE LEFT-HAND SCALE, 
BROKEN CURVES: USE RIGHT-HAND SCALE 
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FIG, 2-2 EQUILIBRIUM COMPOSITION OF AIR AT p/p, = 10-4 
SOLID CURVES: USE LEFT-HAND SCALE 
BROKEN CURVES: USE RIGHT-HAND SCALE 
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FIG, 2-3 EQUILIBRIUM COMPOSITION OF AIR a P/F = 1072 
LID CURVES: USE LEFT-HAND SC. 
BROKEN CURVES; US® RIGHT-HAND pare 
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FIG, 2-4 EQUILIBRIUM COMPOSITION OF AIR AT p/P, =1 
SGLID CURVES: USE LEFT-HAND SCALE, 
BROKEN CURVES: USE RIGHT-HAND SCALE 
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FIG. 2-8 EQUILIBRIUM COMPOSITION OF AIR AT p/P, =] 
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_ Appendix A. Composition and Properties of the Atmosphere” 
(written by A.D. Anderson) 


In this Appendix the main characteristics of the atmosphere are 
reviewed. The principal regions of the eerth's atmosphere are given in 


Table A-1. Inasmuch as the dominant process in the atmosphere is mix- 


_ 


» ing up to at least 90 km, the proportions of the major constitvents are 
“constant up. to this altitude. Hence, the mean molecular weight (mass) 
‘remains the same from 0 to 90 km and the region is called the “homosphere“. 
_ Above 90 km, the molecular weight decreases as the composition changes 
‘with altitude, beceiie of molecular dissociation and diffusion. This 


region is called the “heterosphere”. 


Sgn AE Atlin r te AOI 


Aside from the compositions cf the homosphere and heterosphere, ; 
we have els> summarized the physics of the upper atmosphere with 

| particular attention to variation of properties and composition with the 
sOiar cycle, a subject which has received much recent study. 

Composition of the homosphere (lower atmosphere). The composition : 
of cry air near sea level is presented in Table A-2 (U.S. Standard Atmosphere, | 
1962). Table A-3 gives property values for the homosphere as taken from . <a 
the U.S. Standard Atmosphere, 1962 (1962). The molecular weight from : 
0 to 90 km is 28.964. The U.S. Standard Atmosphere, 1962, is a middle- 
latitude (approximately 45°) year-round mean over the range of solar 
activity between sunspot minima and sunspot maxima. Seasonal and 


latitudinal] variations of the lower atmosphere properties can be taken into 


Lite be 


account, if desired, by using values from supplemental atmospheres derived 
by Cole and Kantor (1963). 


Some of this material has appeared in Space Materials Handbook, Ch. 4, 
published by Addison-Wesley Pub. Co., 1965, edited by C.G. Goetzel, 
J.B. Rittenhouse and J.B. Singletary. 
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TABLE A~1. 


Principal regions of the earth's atmosphere 


Atmospheric 
region 


Sub-region 


Approximate 
altitude 
range (km) 


Homosphere 


Heterosphere 


Troposphere 


Stratospherw 


Mesosphere 


Thermosphere 


Exosphere 


Characteristic features 


0-12 


12-50 


50-90 


90-550 


550-60,000 
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Mean molecular weight 
constant; heat transfer by 
convection 


Constant molecular weight; 
increasing temperature, 
region strongly heated by 
both earth infrared and 
solar ultraviolet radiation 


Constant molecular weight; 
decreasing temperature. 
Mixing processes dominant 
throughout homosphere 


Frequent particie collisions; 
diffusion process dominant 


Collisions rare; temperature 
constant with altitude to 
about 8500 km; diffusion 
process dominant; heat 
transfer by conduction. 


TABLE A-2. 


near sea level from U. S, Standard Atmosphere, 1962. 


Normal composition of clean, dry atmospheric air 


Constitutent gas 


Content 
(per cent by volume) 


Molecular weight 
(c}? = 12,0000) 


Molecular nitrogen 78.08) 26.0134 
Moleculer oxygen 20.9476 31.9988 
Argon 0.93b 39.948 
Carbon dioxide 0.0314 4.00995 
Neon 0.001818 20.163 
Heliym 0.00052h h.026 
Krypton 0.000114 83.60 
Xenon 0.0000087 131.30 
Molecular hydrogen 0.00005 2.0159k 
Me thane 0.0002 16.04303 
Nitrous oxide 0.00005 4.0128 
Sulfur dioxide 0 to 0.0001 64.0628 
Ozone Summer: 0 to 0.000007 47.9982 
Winter: O to 0.000002 
Nitrogen dioxide 0 to 0.000002 46.0055 
Ammonia 0 to trace 17.03061 
Carbon monoxide 0 to trace 26.01055 
Iodine 0 to 0.000001 253.8088 
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Minor constituents in the homosphere. Despite the general constancy 
of the proportions of the major constituents in the homosphere, photochemical 
and collision mechanisms result in major changes with altitude in the concentrations 
of the minor constituents. In spite of their small concentraticas, some of 
these minor constituents can be very important. For example, although water 
vapor usually comprises less than 3 per cent of the gases even with moist 
conditions at sea level, it absorbs nearly six times as much solar radiant 
energy as do all the other gases combined. Furthermore, it accounts for 
nearly all the gaseous absorption of the terrestrial infrared radiation. 

Atomic oxygen is present as a minor constituent of the mesosphere, as 
a result of photodissociation. Detaction of the sodium D lines in the night 
airglow and twilight flash demonstrates the presence of sodium in the mesosphere. 
Carhon dioxide and the oxides of nitrogen have also been detected in the mesosphere 
by mass spectrometers of rocket flights. The presence of water vapor in the 
stratosphere and mesosphere is inferred from the detection of the hydroxvl 
radical in the night airglow. 

Carbon dioxide is also an important absorber of infrared radiation. 
Many measurements have been made to determine the co, content of the 
atmosphere. Bray (1959) has weighted the various measvrements and gives 
a median value of 320 ppm. Most of the carbon dioxide measurements have 
been taken at sea level. Glueckauf (1944) reports values that vary from 250 ppm 
to 300 ppm for a series of balloon flights over England. He indicates that 
samples taken in the stratosphere were not different from those taken at 
ground level. Measurements indicate that the minimum concentrations of co, 
are approximately 150 ppm and are found in polar air; however, the polar air veiges 
vary widely, The maximum concentrations have been measured Off the west coast of 
Africa and may rise to 700 ppm locally. The co, concentration is greater in 
urban than in rural areas, continental air masses show higher concentrations than 


air masses over the Oceans, and night alr us ially contains more co, than daytime air. 
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Water vapor is the most variable constituent in the atmosphere and is 
pezhaps the most difficult to measure at small concentrations and at low temperat 
Measurements indicate that the mixing ratio of water vapor in the atmosphere 
decreases rapidly from ground level to the tropopause. In the stretoaphere, 
the vapor content changes more siowly and usually decreases with altitude. Sinc 
the maximum concentration of water vapor is tempsrature dependent, its normal 
altitude profile is in accord with the temperature profile in the atmosphere. 
Gutnick (1962) has derived a water vapor profile for temperate !atitudes, based 
on the best water vapor measurements to date. The genera! features of this 
model are that the mixing retio decreases from about 6,150 ppm (mg/kg) at the 
surface to 9ppm at 16 km, then increases slighti; with altitude. However, the 
water vapor increase above 16 km is now considered dubious. Bandeen et al. (19 
infer the relative planetary distribution of atmospheric moisture from measurement. 
within the 6.3 micron band of water vapor made from Tiros meteorological 
satellites. 

Oxone, produced by the action of solar ultraviolet radiation on molecular 
oxygen, is found between 15 and 40 km, witn peak values between 22 and 35 km. 
The total amount of ozone in temperate latitudes is about 2.7 mm on the average 
for the whole year; that is, if the total amount of ozone in the atmosphere could 
be concentrated at the surface it would make a layer of this height at STP. 

There are cenaiderable changes in the amount of ozone fron day to day which are 
correlated with weather changes. In particular the relatio1 to pressure changes 
is well indicated. Ozone also appears to have a slight seasonal! and latitudinal 
variation. Figure A-1 shows the measurements for various latitudes reported by 
Miller (1960). The ozone distributions shown in the figure can be considered 
as representative of the atmosphere over the Tropical, Temperate, and Arctic 
reyions of the earth. Atmoapheric absorption by ozone is moat important in the 
ultraviolet region of the spectrum and near 9.6 microns in the infrared; a minor 
absorption band also occurs near 4.8 microns. A good review article on ozone 


has been presented by Gétz (1951). 
316 


ES SY LA ACE LIE: INS A ATL PT Ee eR SESE a 


Altitude (km) 


50 


30 


ALCTIC 


TROPICAL — KODAIKANAL (10° 14! N) 


ARCTIC ~ TROMSO (69° 40' N) 


5 10 
OZONE CONCENTRATION (Atms-cm x 1072 km) 


FIG. A-1 OZONE DISTRIBUTION 
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Gompositioa of the heterosphere (upper atmospnern:. The neutral 
atmosphere above 100 km is composed almvst whclly of molecular nitrogen 
and oxygen and atomic nitrogen, oxygen, halium, and hydrogen; the 
relative concentrations of these constituents depend strongly on altitude 
and temperature. The composition of the upper atmosphere can be explained, 
ai least in a qualitative sense, by noting the types of photochemical reactions 
that can occur and how the reaction products are affected by mixing and 
diffusion. In the following paragraphs, the reactions leading to the neutral 
constituents will be discussed first, followed by the fon reactions. 

At about 100 km, the absorption of solar radiation with wavelengths 
shorter than 1850A down to about 1300A leads to the dissociation of oxygen 
molecules into oxygen atoms. For wavelengths shorter than 1026A, the 
oxygen molecule can he ionized; this ionization is normally followed by a 
dissociative recombination producing atomic oxygen. Although the oxygen 
atoms can recombine into molecules, photochemical equilibrium, where the 
rates Of dissociation and recombination are equal at each altitude, does not 
prevail. Instead, vertical transport due to both diffusion and mixing plays 
an important role in determining the atomic and molecular concentrations 
at various levels near 100 km. More oxygen dissociates than recombines 
above 100 km, due to the rapid fall~off with altitude of the recombination 
processes. Below 100 km, collisions occur frequently enough for recombina- 
tion to prevail, hence more oxygen recombines than dissociates. Conse~- 
quent)y, there is a steady flux of molecular oxygen upward and atomic 
oxygen downward through the 100 km level, due to the effects of diffusion 
and mixing. Atomic cxygen is the most important conatituent in the upper 
thermosphere. 

The most active process leading to the dissociation of molecular 


nitrogen is ionization followed by dissociative recombination, producing 
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atomic nitrogen, Atomic nitrogen can react with molecular oxygen to form 
nitric oxide and atomic oxygen. The nitric oxide in turn reacts with atemic 
nitrogen to form molecular nitrogen and atomic oxygen. The effectiveness 

of these reactions, together with the slowness with which molecular nitrogen 
dissociates, causes atmospheric nitrogen to remain predominantly in 
molecular form. 

Photodissociation of water vapor and methane near 8G km constitute 
the principal sources of atomic hydrogen. Owing to the small mass of the 
hydrogen atom compared to other atmospheric constituents, the hydrogen 
concentration does not decrease with altitude as rapidly as do the other 
atmospheric constituents in the altitude region where diffusion proceeds 
rapidly; hence atomic hydrogen becomes an increasingly important atmos- 
Pheric constituent with increasing altitude. However, atomic hydrogen is 
such a minor constituent in the thermosphere that it does not become the 
dominant constituent until 2000 to 5000 km altitude, remaining so until 
about 20,000 km where the hydrogen ion becomes dominant. The source of 
atomic hydrogen near 80 km can be expected to remain essentially constant 
through the sunspot cycle, but the rate of escape, depending on the temperature 
at the base of the exosphere, varies with the sunspot cycle. The escape 
will be relatively rapid when the temperature is high, and the concentration 
of hydrogen will be correspendingly low in the exosphere near sunspot 
maximum. The escape is relatively slow when the temperature is low, so 
that the concentration must be comparatively high near sunspot minimum. 

Nicolet (1961) showed that helium atoms are an important constituent 
in tne lower exosphere. He explains the high densities derived from the 
rate of change of the period of the Echo satellite by the presence of helium. 
Evidently, atomic oxygen, nitrogen or hydrogen cannot explain the siow 
density decrease between 750 and 1500 km. Atomic oxygen is the mest 


important constituent in the upper thermosphere, but atomic helium 
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dominates over atomic oxygen somewhere above 800 to 1000 km. Atomic 
hydrogen dominates over helium somewhere above 2000 to 5000 km. 
Although the total concentration of the neutral particles is at least an order 
of magnitude greater than that of the fons up to about 2000 km altitude, 
the ion concentr tions finally become significant at the higher altitudes. 

The primar.’ ions formed in the E region (85 to 140 km) of the 
ionosphere are N3 5 Oo; , and o* » The Nj ions dissociatively 
recombine very rapidly, and they may also react with oxygen, so that 
the ooncentration of N5 is small. Ion-atom exchange reactions of the 
type 

X +yz=xv4z (Al. 


proceed very rapidly, so that the reactions of o* with Oo, and Ny 
quickly remove the O° ions and produce sither NO* or oO; - The 
dominant ions in the F, region (140 to 200 km) are NO” and OF near 

the lower altitude boundary, with a gradual transition to of as the 
Principal ion at the upper boundary. ot becomes the dominant ion 
because of the rapid decrease with altitude of the neutral molecular 
constituents that otherwise would tend to eliminate O* by reaction (Al.). 
In the F. region (200 to 800 km), the ions present are O* and N* , with 
rod greatly predominant. The helium ion starts to become dominant between 
800 and 1400 km (Bourdeau and Bauer, 1963), depending on atmospheric 
temperature. Likewise, protons become the dominant ion between about 
1400 to 4000 km, Thus, there are two transition regions (from oxygen to 
helium ions and from helium to hydrogen ions) in the upper ionosph sre 
rather than a single transition from oxygen to hydrogen as previously 
believed. The outer portion of the earth's atmosphere (extending from 


about 29,000 km to at least 60,000 km above the earth's surface) consists 
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mainly of protons. The protons predominate in the earth's outer 
atmosphere over the hydrogen atoins because the confining effect of 
the earth's magnetic field makes their rate of escape to space much less 
than for the neutral hydrogen atoms. According to Johnson (1960), the 
solar wind cannot penetrate the geomagnetic field; therefore, the upper 
boundary cf the earth's atmosphere and magnetic field (magnetosphere) 
are the same. 

Ever. under normal solar wind conditions, the magnetosphere does 
not present a simple boundary to the earth's atmosphere. On the day side 
of the earth, a cavity, which is filled with the earth's magnetic field and 
atmosphere, is formed in the solar wind. At the same time, the geomagnetic 
field is deformed by the solar wind. A schematic representation of the 
cavity and deformed geomagnetic field during a moderately strong solar wind 
is shown in Fig. A~2 (Johnson, 1960). The boundary of the magnetic field 
in the direction facing the sun is about 10 earth radii from the center of the 
earth. The widest portion of the magnetic field downwind from the earth 
is more than 12 earth radii in width and may be over 18. The length of 
this tail is presently unknown and undoubtedly depends on the intensity and 
temperature of the particles in the solar wind. 

Physics of the atmosphere. Assuming that the earth's atmosphere is 
a continuous medium consisting of a gas in static equilibrium, the equation 


connecting the pressure P and density po at any altitude h is 


qdP = -gedh, 
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where g is the gravitational acceleration. The variation in g with 


altitude can be taken into account by using the relationship 
ee 2 2 
D9 = 9,R,AR, +h)” , (A3.) 


where Io is the acceleration of gravity at the earth's surface and Ry 
is the earth's radius. Eq.(A2.)is the hydrostatic equation in differential 
form. Further, let m be the mean molecular mass of the gas and n 


denote the number density or concentration; the density p is then 


ep =nm. (A4.) 


Since the terrestrial atmosphere can be approximated by a 


perfect gas, the pressure is given by the equation of state 


P=nkT , (AS.) 


where T is the temperature and k the Boltzmann constant (1.38 x 1076 
ly 


. 


erg deg™ 


Multiplying both sides of Eq.(A4.)by a and dividing the result into 


Eq. (A5.)gives 
Ph g = kT/mg =H , (A6.) 


where H=k T/mg _ is known as the pressure scale height, a parameter 
convenient for atmospheric calculations. 


From Eqs.{A2.) to (A4.)one can obtain 


=. mg 
hp dh. (A7.) 
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Then, integrating Eq.(A7.) gives the barometric law 


h 
T 
n . 0 _{mg 
n, ol ao an 


and Eqs.(A4.) and(A9.) yield 


h 
m T 4 
eles | _{ mg 
& mT ox | [23 a 
oO 


In Eqs.(A8.) to(Al0.),P, n, and p are the pressure, particle 
concentration, and dens:ty, respectively, at altitude h, and P, rf. 


and P, are the corresponding values at an arbitrary reterence level 


at which h is assigned the value zevo. 

Eas .(A2.) to(Al0.)are applicadle in the atmosphere where the escape 
into space of a constituent is not important. Ths hydrostatic Eq.{A2.)is 
based on the concept of LTE, i.e., collisions between particles are 
frequent enough so that there is a complete statistical axchange cof 
particle energy and momentum in the volume under consideration. In the 
thermosphere the particles collide often enough to insure a Maxwellian 


velocity distribution and the existence of a meaningful kinetic temperature. 
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In the exosphere, above about 550 km, collisions are sufficiently reduced 
so that the existence of a Maxwellian distribution, and consequently, the 
validity of applying the hydrostatic relationship, must be examined anew. 
Liouville's theorem can be used to show that the Maxwellian velocity 
distribution present below the base of the exosphere applies equally well 
in the exosphere, provided that the escape of particlan to space is 
negligible. The atmosphere is sorarefied in the exosphere that there is 
little solar radiation absorbed. Also, there are no energy loss mechanisms 
to disturb the Maxwellian velocity distribution. At the temperatures 
present, the distribution is disturbed only by the escape of neutral 
hydrogen from well above the base of the exosphere. It follows that the 
maximum altitude in the atmosphere where the hydrostatic equation is 
applicable can be determined by calculating the altitude in the exosphere 
where the density of the neutral hydrogen starts to be significantly 
affected by the escape of hydrogen atoms. Consequently, the greatest 
error involved in using the hydrostatic assumption will occur during the 
day near sunspot maximum conditions, inasmuch as temperatures are 
highest then. Calculations (Anderson and Francis, 1964) indicate that 
the error will start to become significant above about 5000 km altitude. 
For average sunspot conditions, the altitude is about 8500 km. For 
gunspot minimum conditions, the error will not be significant below 10,000 km, 
Inasmuch as the atmosphere consists of a mixture of gases subject 
to a force field, the equilibrium distribution of its constituents can be 
expected to shov' some degree of diffusive separation. Although turbulent 
mixing below 160 km does not allow the development of diffusive equilibrium 


distributions, above this altitude there is experimental evidence that 
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diffusive equilibrium does exist. In the geopotential field, diffusive 
equilibrium gives concentration distributions for neutral particles that 

vary exponentially with the geopotential, with a more rapid decrease in 
concentration with increasing geopotential for the heavier constituents 

in the atmosphere than for the lighter. The result is a static distribution 

of gas constituents under the action of the external force field, which may 

be expressed by the barometric law Eq. (A8.) that applies for each atmospheric 
conatituent independently of ths others. 

The concept that the atmosphere extending above 400 km over a given 
location on the earth's surface is isothermal in the sense that the temperature 
does not vary with altitude is now well established. In the tenuous gas 
of the upper thermosphere the thermal conductivity is independent of the 
pressure while the heat capacity varies linearly with density. Consequently, 
the conductivity is very large compared to the heat capacity. Above 400 km, 
the absorption of energy is neoiigible and the relatively high heat conduc- 
tivity eliminates temperature differences; hence the kinetic temperature ia 
nearly constant with altitude for many thousands of kilometers. The 
kinetic temperature can only be determined for a gas with a Maxwellian 
velocity distribution. As pointed out under the discussion of the accuracy 
of the hydrostatic assumption, the Maxwellian velocity distribution applies 
in the exosphere provided thet the escape of particles to space is negligible. 
For hydrogen, the escape of atoms is comparatively rapid so that the 
velocity distribution in the upper exosphere beyond 10,000 km is not 
Maxwellian. Therefore, the hydrogen atoms in the upper exosphere hava 


a non-Maxwellian distribution that bacomes more proncunced with altitude. 
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Under these circumstances the concept of kinetic temperature ia not 
entirely applicable, although an effective temperature can be defined 
by considering the average energy of the hydrogen atoms (Sharp, 1962). | 
Variation of upper-at re properties and osition wit : 
the solar cycle. A principal problem connected with atmospheric structure 
is to calculate accurate values for all the primary properties in the hetero- 
sphere (above 90 km) as a function of time, location, and date by taking 
into account the relevant processes. The primary atmospheric properties 
are temperature, pressure, density, and mean molecular mass (or 
composition), To account for their variations, the following factors 
must be considered: (a) time (hour, day, sun-rotation period season, 
year, sunspot cycle); (b) location (altitude, latitude, 'ongitude); 
(c) solar characteristics (ultraviolet radiation, solar plasma and associated 
magnetic storms); and (d) particle processes (conduction, diffusion, 
mass transport, photoicnization, dissociation, recombination, particle 
escape into space). The problem of describing upper atmospheric 
behavior is difficult because many of the above elements are interrelated. 
The magnitudes of the atmospheric properties in the upper atmosphere 
are derived from measurements made from satellites, rockets, meteor 
Gbwswationa sky emissions, and the propagation of sound and radio 
waves. The data are sparse and contain much uncertainty above 200 km. 
The density (drag) data resulting from tracking satellites are the most 
precise, and by far, the most numerous. Study of the orbital decay data 
has clearly established that two major systematic density variations 
occur: (1) a solar activity effect in which variations in atmospheric 


heating and density take place above 90 km due to variations in solar 
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ultraviolet radiation; and (2) a diurna! (time-of-day) effect, tn which © 
the solar heating results in the atmosphere bulging toward the sun, 
producing relatively large density increasas at altitudes above 300 km 
in the sunlit region of the earth. At 800 km, due to effect (1), the 
density can be 40 times as high during solar maximum conditions than 
during solar minimum, and due to effect (2), 15 times as high during the 
day than during the night. The combination of both (1) and (2) can 
result in densities 500 times as high at 800 km during solar maximum (day) 
than during solar minimum (night). This extreme variability in density 
also applies to pressure. Upper atmosphere density variations with 
latitude and season are much smaller than the two primary effects above. 
In addition, there is a conspicuous drag effect on satellites associated 
with major magnetic storms and the rasulting heating is such that the 
relative effect on density increases with altitude. At present, it is not 
possible to make more than a rough estimate of the effect of magnetic 
heating on the density; fortunately, major magnetic storms are fairly 
rare and the associated heating is transient. 

Upper-atmosphere measurements are not made on a regular enough 
basis, either in kind, time or space, to allow them to be used alone, 
without the aid of a model, to represent atmospheric conditions. The 
approach followed by most models used to derive atmospheric properties 
ig tc assume altitude profiles for some of them in order to calculate the 
remainder. Almost all the models deal with data referring only to density, 
pressure and/or temperature. The altitude variation of the mean molecular 


mass is introduced somewhat arbitrarily and, therefore, a physically 
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consistent vertical distribution of the composition cannot be obtained. 

A new method has been devised for computing atmospheric properties . 
(Anderson and Francis, 1964). In this model, no major assumptions are 
made regarding the property profiles. Instead, the primary properties 

are calculated by starting with an empirical density profile from a 
density model as the chief input, assuming diffusive equilibrium 
conditions above 110 km and isothermal conditions with altitude above 
400 km. The density profile, used as starting input for this new property 
model, is represented in a previous empirical mode! (Anderson, 1964) as 
a function of local time and solar activity from 200 to 800 km. This 
density model can be used to calculate the marked variations in density 
that occur from day to day by taking corrected values of the 10.7-cm 
solar flux (called S') as an index of the sun's extreme ultraviolet 
radiation. The remaining properties are derived by using the hydrostatic 
equation and equation of state, that relate density, pressure, temperature, 
and mean molecular mass as a function of altitude (Eqs.(A8.)to(Al0.), 
together with boundary values based on measurements. 

The equations embodying the model are programmed for au IBM-7090 

computer, The results of the computations made from the model are 

given in two tables that present the neutral atmospheric properties and 
composition vs. altitude from 100 to 10,000 km. Tables A-4 and A-5 
exhibit the properties and number densities of the various constituents 
vs. altitude near sunspot maximum and near sunspot minimum, respectively. 
The S' value used for Table A-4 (S' = 250) is an average value for 


September, 1959, the month immediately following the Argus high-altitude 
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nuclear detonation. The S' value used for Table A-5 (S' = 44) is an 
average value for July and Augusi, 1962 (S* = 42 for July: S' = 46 for 
August), the two months following the Starfish high-altitude shot. Both 
tables are for t = 21 hr.; the density for this time has been found to 
closely approximate the diurnally averaged density or the sum of the 


densities for every hour of the day divided by 24. Calculations of 


density from the model are based on neutral particles only. The ions do 


not contribute significantly to the total (ambient) density below 2000 km. 


Above this altitude, the contribution of H* should be taken into account 
in any calculation of the ambient density. The neutral density is 
corrected for the departure from a Maxwellian velocity distribution caused 
by the escape of atomic hydrogen from the earth's exosphere. The 
temperature u: ad to calculate the total pressure in the exosphere from 

A-5 is the kinetic temperature that is assumed to be locally isothermal 
vertically to 10,000 km. Hence, the computed pressure will be somewhat 
higher than the ambient pressure (a function of the effective temperature) 


of the neutral constituents above about 5000 km. 
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